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ABSTRACT 


The  solution  of  the  continuous  waveguide  transition  problem  can  be  obtained  by 
discretizing  the  boundary  and  applying  mode-matching  or  by  using  a  system  of  oroinary 
differential  equations.  Both  approaches  involve  approximate  representations  of  the 
boundary.  When  using  the  differential  equation  approach,  it  was  found  necessary  to 
consider  tlie  transition  as  several  sections  in  series  in  order  to  avoid  numerical  instabilities. 
When  this  is  done,  one  may  cascade  using  a  generalized  scattering  matrix  approach  or  a 
generalized  ABCO  matrix  method.  Results  are  shown  comparing  the  accuracy  of  the 
boundary  discretization  approach  and  the  differential  equation  approach  for  the  Marie 
transducer  and  for  linear  transitions  of  various  lengths  in  rectangular  waveguide;  experi¬ 
mental  results  are  also  given  for  the  Marie  transducer. 
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1.  INTRODUCTION 


Numerical  approaches  to  the  continuous  waveguide  uansition  analysis  and  design 
problems  have  been  studied  for  the  past  twenty  or  so  yeais.*'^  Recently,  the  author  of 
this  paper  used  a  novel  moment  method  technique  to  solve  a  system  of  ordinary 
differential  equations  describing  such  a  transition/-^  this  paper  is  basically  an 
extension  of  the  work  described  in  References  4  and  S.  The  paper  compares  the 
differential  equation  method  with  an  approach  that  approximates  the  transition  as  a 
series  of  discrete  steps  and  subsequently  ^plies  mode-matching,  applies  these 
techniques  to  linear  rectangular  transitions,  discusses  the  validity  of  the  differential 
equations  and  assoc  iated  boundary  conditions,  and  reports  experimental  and  numerical 
results  for  the  Marie  rectangular  TEio  to  circular  TEoi  mode  transducer.  Also  included 
are  new  results  regarding  the  sensitivity  of  the  computed  solution  to  the  choice  of 
cascading  formulation  when  (as  is  a  fairly  common  practice)  a  transition  is  analyzed 
as  several  sections  in  series  in  order  to  avoid  numerical  instabilities.^-^-^ 


2.  METHODS  OF  SOLUTION 


The  system  of  differential  equations  used  in  this  paper  is  based  on  the  assumption 
that  the  transverse  portion  of  the  electromagnetic  field  in  a  transition  can  be  written  as 
a  sum  of  the  uniform  waveguide  mode  functions  conesponding  to  the  local  transition 
cross-section 


E,ix,y,z)  =  '^V„(,z)e„ix,y,z)  (1) 

•• 

H,(x,y,2)  =  ]^/„(z)*„(x,y,z)  (2) 

where  z  denotes  the  direcdon  of  propagaUon,  and  the  parameters  Vn(z)  and  I«(z)  are 
known,  respecdvely,  as  the  equivalent  voltages  and  the  equivalent  currents.  A  system 
of  otxiinary  differential  equations  that  specifies  V,„(z)  and  Lfz)  was  given  by  Reiter  in 
19S9,‘  and  methods  of  soluuon  are  discussed  at  some  length  in  References  3,4,  and  S. 
These  equauons  are 


=  +  (3) 

=  — T  I  (4) 

**  H*! 

Both  TE  and  TM  modes  are  included  in  Equauons  3  and  4.  The  variable  0m  denotes 
the  wave  number  of  the  with  mode,  the  variable  Zm  denotes  the  wave  impedance  of  the 
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mth  mode,  and  the  “transfer  coefficients”  Tm  (which  describe  coupling  between  the 
two  inodes  m  and  n)  are  given  by* 

(5) 

SU) 

The  integration  in  Equation  S  is  over  the  local  waveguide  cross-section  S(z).  If  the 
system  of  Equations  in  3  and  4  is  truncated  and  then  solved  for  the  portion  of  the 
transition  between  z  =  z,  and  z  =  Z{«i,  one  may  specify  numerical  values  for  the 
associated  generalized  ABCD  matrix 


■V'(Z.>i)' 

Bi' 

■V(Z.)‘ 

./(Z.>i). 

A. 

./(Z.), 

where  the  dimension  of  the  square  submatrices  A,,  B,.  C„  and  D,  is  equal  to  the  number 
of  modes  (propagating  and  evanescent)  included  in  the  solution.  Cascading  the  results 
for  two  adjacent  portions  of  the  transition  may  be  accomplished  by  simple  multiplica¬ 
tion  of  the  two  ABCD  matrices,*'^  or  by  translating  these  two  matrices  into  gener^ized 
scattering  matrices  and  then  linking  them  together.^  One  goal  of  this  paper  is  to 
compare  the  efficacy  of  these  two  cascading  approaches.  The  ABCD  matrix  approach 
is  computationally  less  intensive  and,  therefore,  initially  more  attractive  than  the 
scattering  matrix  approach. 

One  well-known  alternative  to  solving  the  differential  equations  is  to  discretize  the 
boundary  and  to  apply  mode-matching.  The  boundary  discretization  geometry  used  in 
this  paper  is  depict^  in  Figure  1 .  The  two  uniform  waveguides  in  Figure  1  have  cross- 
sections  equivalent  to  the  uansition  cross-sections  at  z  and  z  +  AL.  As  shown  in 
Appendix  A,  the  difference  equation  for  the  equivalent  voltages  and  currents  is 


V^i2  +  AL)  = 


- ■' 2, (z)sin^^i^ j  jje^(z)  e„{z  +  AL)dxdy 


+yZ„(z-^AL)sin(^^j|; 


jje^(z)e„{z  +  M,)dxdy 
S 


^.(z) 
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-/Z,(2  -h  A£,)sin^^^  j 


(  \ 

jje^(^)-^m(^  +  ^)dxdy  I,{z) 

.S 


(7) 


where  the  area  of  integration  S  is  equal  to  the  intersection  of  S(z)  and(z + AL).  Rewriting 
Equation  7  as  a  sum  of  powers  of  ^  and  retaining  only  terms  of  order  (MJf  or  lower, 
one  obtains 


v^(2+m=v^(2) 


-jALZ^i2)P^i2)U(2) 


jj<U)eMdxcfy 


V,{z) 


4 


(8) 


(AL) 


2  - 


•>1 


// 


eZiz)  ■e,{2)dxdy 


V.Cz) 


<1-1 ^ 


+  Z^{2) 


jle'mU)%iz)<ixdy 

VS 


Uz) 


where  the  primes  refer  to  differentiation  with  respect  to  z.  A  similar  equation  exists  for 
In-  Gearly,  in  the  limit  of  vanishing  AL.  Equation  8  reduces  to  Equation  3.  Therefore, 
one  might  expect  that  the  boundary  discretization  method  should  converge  to  the  same 
solution  as  does  the  diffoential  equation  method.  Convergence  rates  can  be  affected 
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by  the  (ALf  terms,  however,  and  these  terms  do.  in  fact,  differ  between  Equation  8  and 
some  of  the  techniques  fpr  solving  Equations  3  and  4  (e.g.,  the  Runge-Kutta  method, 
discussed  in  Appendix  B.  Similar  difference  equations  can  be  derived  when  traveling 
waves  areconsidered  instead  of  the  equivalent  voltages  and  cuirentsofEquations  1  and 
2  (further  details  are  given  in  Appendix  Q.  The  mode-matching  ^proach  used  in  this 
paper’  is  a  modification  of  the  method  described  by  Carin,  Webb,  and  Weinreb.*  The 
scattering  matrix  cascading  technique  described  by  Chu  and  Itoh^  is  used  to  link 
together  the  large  number  of  junctions  used  to  approximate  a  transition. 

Although  Equations  3  and  8  indicate  that  the  differential  equation  method  and  the 
boundary  discretization  method  should  converge  to  a  common  solution,  the  two 
methods  are  based  on  somewhat  different  representations  of  the  transition  boundary. 
The  differential  equation  approach  assumes  a  perfectly  smooth  boundary,  while  the 
mode-matching  approach  uses  a  stepped  boundary.  It  may  not  be  correct  to  say  that  the 
stepped  boundary  is  a  valid  approximation  of  the  smooth  boundary.  This  is  well 
illustrated  if  one  attempts  to  approximate  the  diagonal  of  a  unit  square  as  a  series  of 
steps;  the  approximate  boundary  will  have  a  length  equal  to  two,  independent  of  the 
number  of  steps,  whereas  the  actual  diagonal  has  a  length  equal  to  .  Nonetheless, 
in  spite  of  these  two  rather  different  boundary  representations.  Equations  3  and  8 
appear  to  indicate  that  the  two  methods  are  equivalent.  This  can  perhaps  be  explained 
as  follows.  Most  standard  methods  for  solving  Equation  3  can  be  written  as  difference 
equations  similar  to  Equation  8.  Any  difference  equation  approach  for  solving 
Equation  3  will  include  terms  of  (AL)^  and  higher,  and,  for  nonzero  AL,  these  higher 
order  terms  will  represent  a  distortion  of  the  smooth  system  of  Equation  3.  This 
distortion  vanishes  as  AL  vanishes,  resulting  in  the  equivalence  of  Equations  3  and  8. 
Yet  another  interesting  issue  is  whether  or  not  Equations  1  and  2  are  appropiate  for 
nonuniform  waveguides.  This  will  be  discussed  in  Section  4.0. 


Z  +  AU2 


Discretized  boundary 


Actual  boundary 


Figure  1.  Boundary  discretization  geometry. 
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3.0  MARIE  TRANSDUCER  ANALYSIS 


In  this  section,  the  numerical  techniques  described  previously  are  aqrplied  to  an  X- 
band  model  of  the  Marie  rectangular  TEio  to  circular  TEoi  mode  transducer,'*’  and  the 
resultant  data  are  compared  with  experimental  results.  This  device,  which  has  been 
informaUy  described  as  “the  best  transition  to  TEoi.""  is  shown  in  Figure  2,  and  the 
desired  electric  field  lines  at  various  locations  are  shown  in  Figure  3 .  (The  circular  taper 
portion  in  Fig.  2  is  not  included  in  some  models.) 

A  detailed  study  of  the  Marie  transducer  was  performed  by  Saad,  Davies,  and 
Davies  in  the  1970s.'‘  ^  One  dinerence  between  their  work  and  ours  is  that  they 
numerically  implemented  Solymar’s  small  coupling  procedure,'^  while  we  use  the 
methods  described  in  References  4  and  5.  A  proequisite  for  analyzing  a  continuous 
waveguide  transition  is  to  obtain  the  mode  functions  in  Equations  1  and  2.  These  modal 
solutions  and  the  associated  eigenvalues  are  essential  in  determining  the  coefficients 
in  Reiter’s  equations.'*'*  As  noted  in  Reference  12,  these  tasks  must  be  accomplished 
numerically  when  one  considers  an  irregularly  shaped  device  such  as  the  Marie 
transducer.  The  publicly  available  software  we  used  is  described  in  References  14  and 
IS,  and  further  details  regarding  our  particular  implementation  can  be  found  in 
Reference  7.  In  the  circular  taper  portion  of  the  transition,  the  coefficients  for  the 
differential  equations  may  be  determined  analytically.'* 

The  generation  of  numerical  data  describing  the  transmission  and  reflection 
properties  of  this  device  will  now  be  described.  These  data  will  be  compared  with 
measurements  performed  on  two  Marie  transducers  connected  to  each  other  by  a  35- 
inch  long  metallic  circular  waveguide,  i.e.,  launching  and  detecting  at  the  two 
rectangular  ports.  Each  Marie  transducer  was  46.35  inches  long  with  circular  taper 


Figure  2.  The  Marie  transducer. 
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[] 

HDD 

Figur*  3.  Desired  electric  field  lines  for  the  Marie  transducer. 

portions,  each  11  inches  long.  The  circular  ends  of  the  Marie  transducer  were  1.18 
inches  in  radius,  and  the  rectangular  ends  were  standard  WR-90  waveguide.  The 
numerical  data  were  acquired  using  the  following  four  schemes: 

Approach  7 .  In  the  irregularly  shaped  portions  of  the  Marie  transducer,  only  the 
desired  rectangular  TE,, —  circular  TE„  mode  is  included  in  the  expansions  of 
Equations  1  and  2,  and,  in  the  most  narrow  part  of  the  circular  taper,  only  die  desired 
TEg,  mode  was  considered.  Only  a  single  mode  was  considered  because  previous  work 
indicates  that  considering  only  modes  near  or  above  cutoff  is  sufficiently  accurate 
while  it  reduces  computational  complexity.*  All  of  these  "one-mode"  sections  were 
analyzed  using  the  Runge-Kutta  technique  to  solve  the  truncated  system  of  differential 
equations.*  In  a  middle  1  -inch  portion  of  the  taper,  where  the  unwanted  TE,,  mode  goes 
through  cutoff,  both  this  mode  and  the  TE„  mode  were  included  and  the  moment 
method  of  Reference  4  was  used.  The  reason  for  selecting  this  technique  is  that,  as 
observed  in  References  3  and  4,  evanescent  modes  cause  numerical  problems,  and  the 
results  reported  in  Reference  4  seem  to  indicate  that  the  moment  method  technique  has 
better  stability  characteristics  than  do  some  conventional  techniques  (e.g.,  the  Runge- 
Kutta  method  and  an  iterative  integration  technique*-’)-  Finally,  in  the  widest  part  of  the 
taper,  where  both  modes  were  above  cutoff,  the  numerically  less  expensive  Runge- 
Kutta  technique  was  used.  As  in  Reference  4,  cascading  was  accomplished  through  the 
simple  multiplication  of  the  generalized  ABCD  matrices  representing  adjacent  cross- 
sections.  First,  the  4  X  4  TE,, — TE,,  ABCD  matrices  representing  the  circular 
waveguide  and  the  wider  portions  of  the  circular  tapers  of  the  two  Marie  transduco-s 
were  combined  to  produce  a  4  x  4  ABCD  matrix.  The  four  elements  of  this  matrix 
describing  TEg,  to  TEg,  interaction  were  used  to  produce  a  2  x  2  ABCD  matrix,  which 
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was  then  cascaded  with  the  “single-mode”  sections.  This  new  2x2  matrix  was  then 
translated  into  a  2  x  2  scattering  matrix.  (For  further  details,  see  Chapter  4  of  Ref.  7.) 

In  order  to  verify  these  calculations,  measurements  were  performed  using  a 
Wiltron  360  Network  Analyzer  set  to  operate  between  11  and  12  GHz.  These  tests 
iiKluded  a  calibration  procedure  that  enabled  normalization  of  the  test  data  in  the 
presence  of  imperfect  connections,  cables,  etc.,  between  the  unit  under  test  and  the 
network  analyzer  input  and  output  ports.  (To  further  establish  the  validity  of  these  tests, 
all  measurements  were  repeated  several  times,  and  tests  wa%  performed  on  known 
devices.  We  believe  that  the  IS21I  data  are  accurate  to  within  ±0.01  dB  and  IS  1  il  to  within 
±10.0  dB.  (Fw  funner  details,  see  Appendix  D.)  Figure  4  shows  ni  “.rical  and 
experimental  data  for  the  scaaering  parameter  S21  for  the  two  Marie  transducers  and 
the  connecting  circ ular  waveguide.  Figure  5  shows  the  results  for  S 1 1 .  The  experimental 
data  indicate  very  low  but  observable  dissipative  loss. 

Approach2.  The  second  method  of  analysis  differed  firom  Approach  1  in  only  one 
respect:  cascading.  The  individual  4x4  ABCD  matrices  were  converted  into  4X4 
scattering  matrices  before  cascading.  Chu  and  Itoh  have  published  a  method  fex- 
cascading  two  devices  represented  by  generalized  scattering  matrices  and  connected 
by  a  uniform  waveguide.’  In  cascading  two  adjacent  portions  of  the  transition,  this 
method  was  used  with  the  length  of  the  uniform  waveguide  set  equal  to  zero  in  the 
computer  program.  Once  the  “two-mode”  sections  were  cascaded,  the  four  elements 
of  the  generalized  scattering  matrix  representing  TEoi  to  TEoi  interaction  were  used  to 
produce  a  2  x  2  scattering  matrix,  which  was  then  cascaded  with  the  “single-mode” 
sections.  Cascading  such  2x2  matrices  is  a  trivial  task,  and  was  done  according  to  a 
technique  described  in  Reference  17.  Figures6  and  7  show  the  data  generated  using  this 
approach  along  with  the  experimental  data  from  Figures  4  and  5. 


Frequency  (GHz) 

Figur*  4.  Experimental  and  numerical  IS^J  for  system  consisting  of  two  Marie 
transducers  connected  by  a  circular  waveguide.  These  data  were  generated  using  the 
ABCD  matrix  cascading  formulation  (Approach  1  in  text). 
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Figure  5.  Experimental  and  numerical  |S„|  for  system  consisting  of  two  Maris 
transducers  connected  by  a  circular  waveguide.  These  data  were  generated  using  the 
ABCD  matrix  cascading  formulation  (Approach  1  in  text). 


Approach  3.  The  third  approach  differed  from  Approach  2  in  only  one  respect:  the 
middle  part  of  the  circular  taper  was  analyzed  using  boundary  discretization  and  mode¬ 
matching.  Five  modes  were  considered  at  five  junctions  and  the  resultant  scattering 
matrix  was  converted  into  a  two-mode  scattering  matrix  (by  simple  truncation)  before 
cascading  with  the  other  “two-mode”  sections.  Results  obtained  using  this  approach 
were  identical  to  those  of  Approach  2,  (i.e..  Figures  6  and  7). 

Approach  4.  Boundary  discretization  and  mode-matching  were  used  to  analyze 
the  entire  circular  taper.  Nine  modes  were  used  at  each  junction  and  the  tap^  was 
approximated  by  fifty  junctions.  A  2  x  2  scattering  matrix  was  extracted  from  the  results 
and  cascaded  with  the  irregularly  shaped  “single-mode”  sections.  These  results  were 
very  close  to  those  of  Approaches  2  and  3  (i.e..  the  thin  lines  in  Figures  6  and  7). 

Figures  4  through  7  show  that,  at  least  for  this  example,  the  computed  solution  is 
insensitive  to  whether  the  differential  equation  technique  or  the  mode-matching 
technique  is  used.  However,  this  is  not  true  regarding  the  choice  of  cascading  method. 
One  might  conclude  that  the  ABCD  matrix  multiplication  technique  used  in  Approach 
1  provided  a  more  accurate  solution  than  the  scattoing  matrix  cascading  technique 
used  in  the  other  three  approaches,  because  the  numerical  and  experimental  curves  are 
more  closely  coincident  in  Figures  4  and  S  than  in  Figures  6  and  7.  It  should  be  noted, 
however,  that  these  calculations  do  not  take  into  account  conductive  losses;  a  “good” 
solution  would  therefore  show  insertion  losses  slightly  less  than  those  in  the  experi¬ 
ment  If  observable  losses  are  expected,  the  scattering  matrix  cascading  data  of  Figures 
6  and  7  appear  more  credible  than  the  Approach  1  data  of  Figures  4  and  S.  Further,  in 
this  case  the  scattering  matrix  formulation  involved  manipulating  2x2  submatrices 
with  condition  numbers  in  the  range  lOP  to  10^,  while  the  generalized  ABCD  matrix 
formulation  involved  manipulating  4x4  matrices  with  condition  numbers  of  10^.  To 
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Figure  6.  Experimental  and  numerical  IS,,!  for  system  consisting  of  two  Marie 
transducers  connected  by  a  circular  waveguide.  These  data  were  generated  using  the 
scattering  matrix  cascading  formulation  (Thin  line:  Approaches  2  and  3  in  text;  dotted 
line:  Approach  4  in  text.) 


Frequency  (GHz) 


Figure  7.  Experimental  and  numerical  |S„|  for  system  consisting  of  two  Marie 
transducers  connected  by  a  circular  waveguide.  These  data  were  generated  using  the 
scattering  matrix  cascading  formulation  (Thin  line:  Approaches  2  and  3  in  text;  dotted 
line:  Approach  4  in  text.) 
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sum  up,  for  this  particular  example,  the  scattering  matrix  formulation  utilized  smaller 
and  bett£7  conditioned  matrices  and  yielded  physically  more  reasonable  results  than 
does  the  ABCD  matrix  formulation. 

A  comparison  of  computer  run  times  for  the  moment-method  program  and  the 
mode-matching  program  used  in  the  preceding  examples  indicates  no  substantial 
advantage  from  the  viewpoint  of  computational  efficiency.  A  similar  conclusion  is 
reached  when  these  programs  are  used  to  analyze  a  simple  rectangular  taper.^  It  should 
be  noted  that  neither  of  these  two  programs  has  yet  been  optimized  with  respect  to  either 
run  time  or  storage.  The  Runge-Kutta  technique  is  by  far  the  least  expensive  of  the 
methods  considered  here,  and,  as  indicated  in  Reference  4,  it  is  sufflciently  accurate  in 
many  cases. 


4.  CONTINUOUS  VERSUS  DISCRETE  FORMULATIONS 


In  Section  3,  our  particular  implementations  of  the  differential  equation  approach 
and  the  boundary  discretization  approach  yielded  virtually  identical  results  when 
applied  to  the  circular  taper  potion  of  a  Marie  transducer.  This  taper  varies  slowly  with 
respect  to  wavelength.  It  also  has  been  reported  that  the  two  approaches  yield  identical 
data  for  a  slowly  varying  rectangular  taper^  and  several  extremely  slowly  varying 
circular  tapers.^  In  this  section,  we  show  that  the  routines  are  also  consistent  for  rapidly 
varying  transitions.  Our  treatment  starts  with  the  rectangular  taper  of  Figure  8  with  the 
taper  length  L  as  a  parameter.  This  taper  is  analyzed  for  extremely  short  lengths  L  using 
the  previously  described  techniques.  The  L  -  0  case,  which  is  simply  an  abrupt 
discontinuity,  is  analyzed  using  a  single-step  mode-matching  routine.  (The  differential 
equation  technique  clearly  will  not  work  forL  =  0.)  The  idea  is  to  see  whether  or  not, 
for  diminishing  lengths  L,  the  waveguide  transition  data  approaches  the  abrupt 
discontinuity  results  and  to  compare  the  differential  equation  data  to  the  boundary 
discretization  data  for  these  nonzero  lengths  L. 

Figure  8  shows  a  transition  between  two  rectangular  waveguides  with  different 
heights.  For  the  dimensions  given,  we  shall  consider  the  band  between  1 .0  and  1 .8  GHz, 
where  only  the  TEio  mode  is  capable  of  propagation.  The  abrupt  (L  =  0)  case  was 
analyzed  using  two  methods,  namely,  a  nine-mode  mode-matching  solution  and  a 


Figure  8.  Steep  rectangular  waveguide  taper. 
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simple  procedure  used  data  from  Figure  S.26-3  in  Marcuvitz’s  book.*'  Figure  9  shows 
the  dominant  mode  reflection  coefficient  Su  obtained  by  these  two  methods  and  that 
the  two  sets  of  data  are  coincident.  S  ubsequently ,  results  were  generated  for  the  lengths 
L  =  O.OSXo,  0.025X0, 0.OOlXo,  and  lO'^Xo,  where  Xo  is  the  firee-space  wavelength  at  1 .0 
GHz.  These  results  were  generated  taking  into  account  the  nine  lowest  modes  and  by 
using  two  analytical  techniques.  First,  boundary  discretization  and  mode-matching 
were  applied  with  100  junctions  betweenz = Oand  z =L.  Second,  the  Galerkin’s  method 
version  of  the  moment-method  technique'  was  used  with  two  cascaded  sections 
between  z = 0  and  z  s  L,  and  with  five  triangle  weighting  functions  per  section.  For  each 
length,  both  sets  of  curves  were  indistinguishable;  these  curves  are  shown  in  Figure  10. 
The  results  for  L  =  10~%i  are  very  close  to  the  L = 0  results  shown  in  Figure  9  generated 
by  using  a  nine-mode,  single-step,  mode-matching  routine. 

For  the  short  transition  example  just  described,  the  differential  equation  technique 
and  the  boundary  discretization  technique  are  in  agreement,  just  as  they  were  for 


Frequency  {GHz) 

Figure  9.  Computed  reflection  coefficient  for  the  steep  taper  of  Figure  7  for  the  abrupt 
discontinuity  (L  >  0)  case.  The  solid  line  was  generated  using  mode-matching  with  nine 
modes.  The  dots  were  generated  using  Fig.  5.26-3  from  Reference  18. 
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gradual  t^rs.  Indeed,  consistency  between  these  two  methods  is  what  one  would 
expect  after  examining  Equations  3  and  8.  This  agreement  between  the  two  techniques 
might  be  theoretically  important  This  is  because,  as  observed  in  Refnence  19,  the 
magnetic  field  expansion  (Equation  2)  cannot  satisfy  the  transition  wall  boundary 
condition  exactly.  At  the  wall,  the  magnetic  field  must  satisfy  (see  Fig.  1 1): 

Hin  =  0  (9) 

The  individual  terms  of  the  expansion  satisfy: 

11  =  0  (10) 

(The  boundary  condition  for  the  electric  field  is  satisfied  by  the  individual  terms  of 
Equation  1.)  If  the  transition  is  a  slowly  varying  (i.e.,  nearly  uniform)  one.  Equations 


Figura  1 0.  Computed  reflection  coefficient  for  the  steep  taper  of  Figure  7  for  the  cases 
L  -  O.OSXg  (dots),  L  «  0.025\  (dash),  L  -  0.001  \  (chaindot),  and  L  - 1 0-*A,j  (chaindot), 
where  is  the  free-space  wavelength  at  1.0  GHz.  Each  of  these  four  sets  of  data  was 
generated  using  the  differential  equation  method  and  the  boundary  discretization 
method  (nine  modes).  The  two  methods  yielded  identical  results. 


12 


The  Johns  Hopkins  Uniweisiqt 

Applied  Physics  Laboratory 
Laurel,  Maryland  2072S6099 


Figure  11.  Waveguide  transition  geometry  for  electromagnetic  boundary  conditions. 

9  and  10  are  approximately  the  same,  but  for  rapidly  varying  tapers,  the  unit  vectors  m 
and  n  become  less  coincident,  the  approximation  Imeaks  down,  and  one  might  question 
the  validity  of  the  differential  equation  method.  For  the  steep  taper  in  Figure  8, 
however,  the  differential  equation  method  yields  the  same  results  as  the  boundary 
discretization  technique,  and  these  two  idential  sets  of  results  converge  to  the  abrupt 
discontinuity  solution  for  diminishing  valuesof  L.  Moreover,  consistency  among  these 
same  techniques  when  applied  to  a  steep  circular  waveguide  taper  has  also  been 
repcHted  (see  Appendix  E).  Thus,  the  differential  equation  ftMmulation  appears  to  be 
suitable  even  for  rapidly  varying  ugiers.  One  possible  explanation  for  this  may  be  fourxl 
in  an  idea  advanced  by  Unger,'’  namely,  that  the  expansion  Equation  2  is  valid,  that  it 
converges  to  the  correct  field  values  inside  but  not  on  the  boundary,  and  that  this 
expansion  therefore  suffers  a  discontinuity  on  the  boundary.  Recently,  it  has  been 
argued  that  the  mathematics  of  the  derivation  of  Reiter’s  equations  are  fully  consistent 
with  the  existence  of  an  expansion  with  these  attributes  (Ref.  7). 
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5.  CONCLUSION 


Several  waveguide  transitions  have  been  analyzed  using  two  different  methods, 
discretizing  the  transition  boundary  and  using  mode-matching,  and  by  solving  a  system 
of  ordinary  differential  equations.  For  gradual  transitions,  the  two  methods  lead  to 
identical  results  and  use  comparable  amounts  of  computer  time  (assuming  the  same 
number  of  modes  is  used  in  both  cases).  For  extremely  steep  transitions  (excluding  the 
step  transition),  our  results  still  indicate  agreement  between  the  two  methods.  Our 
results  also  indicate  that  the  scattering  matrix  cascading  method  (previously  used  by 
Flugel  and  Kuhn^)  leads  to  results  that  are  more  credible  than  those  obtained  when 
cascading  is  performed  by  multiplying  generalized  ABCD  matrices.  The  differential 
equation  technique  may  be  more  generally  applicable  for  continuous  transitions  in  that 
artificial  discontinuities  are  not  introduced,  and  an  arbitrary  number  (including  zero) 
of  evanescent  modes  can  be  introduced. 
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APPENDIX  A. 

DERIVATION  OF  EQUATION  7 


In  Fig.  A-1,  the  electromagnetic  fields  in  the  z-plahe  are  given  by: 

E(z)  =  j^V;,(2)e,{2)  (A-1) 

•=1 

H(2)=£/.(2)h.(2)  (A-2) 

*=1 

Similarly,  the  fields  in  the  (2  +  AL)  plane  are  given  by: 

E(2  +  AZ.)  =  ]^V,,(2  +  AL)e,,(2  +  A£.)  (A-3) 


H(2  +  AL)  =  5^/.(2  +  AL)Ii,(2  +  A£,)  (A-4) 

Msl 

By  applying  the  ABCD  matrix  for  a  uniform  waveguide,  the  equivalent  voltages  and 
currents  at  2  -i-  AL/2  are  seen  to  be  related  to  Equations  A-1  and  A-2  by: 


t]  "  (z)-  (z)  (A-5) 

z  +  AL/2 


FIgur*  A-1.  Boundary  discretization  geometry. 
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The  equivalent  voltages  and  currents  on  the  left  sides  of  Equations  A-S  and  A-6  are  the 
coefficients  of  the  mode  functions  e„(z)  and  h;,(z).  When  expanding  the  fields  at  z-i-  AL/ 
2  in  terms  of  the  mode  functions  +  AL)  and  hm(z  +  AL),  we  shall  denote  the 
corresponding  coefficients  (z  +  AL/2)  and  /a,(z  +  AL/2).  These  new  variables  are 
related  to  those  of  Equations  A-3  and  A-4  by: 

V„(z  +  AL)  =  cos(^^  jv’.f^z +  jZ,(z  +  AL)sm[^^]/.(z  + -y] 

(A-7) 


/^(z  +  AL)  = 


-J 


Z(z  +  AL) 


sm 


2  +  - 


(A-8) 


Enforcing  the  field  continuity  condition  at  z  AL/2,  one  obtains; 


X  +  ■y]h,(z)  =  X/.(z  +  ^ 


Furthermore,  utilizing  the  well  known  orthonormality  relationships  among  the  mode 
functions. 


The  remainder  of  the  derivation  is  straightforward  but  tedious.  First,  one  substitutes 
Equations  A-S  and  A-6  into  Equations  A- 11  and  A- 12,  thereby  obtaining  new 
expressions  for  V„{z  +  AL/2)  and  I„{z  +  AL/2).  These  new  expressions  are  substituted 
into  Equation  A-7,  immediately  yielding  Equation  7  6rom  the  main  text. 


A-2 
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APPENDIX  B, 

APPLICATION  OF  THE  RUNGE-KOTTA  METHOD  TO  EQUATION  3 


Equations  3  and  4  from  the  main  text  may  be  written  symbolically  as 


where: 


and 


^  =  a2)y 

dz 


V'2(2) 

/,(0 

W 


(B-1) 


(B-2) 


(B-3) 


The  Runge-Kutta  method  may  be  summarized  as  follows.  We  deHne  the  four 
vectors: 


X,  =ALC{z)y(z) 


(B-4) 


^2 


(B-5) 
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(B-6) 

=  ALcj^z  -t-^j(y(2)  +  ATj)  (B-7) 

The  estitnate  for  y(z  -i-  AL)  is: 

y{i^  M.)  =  y{z)+-(K^-^zK-i  +  iK-i^  K^)  (B-8) 

6 

Expanding  the  matrix  C  in  the  form  of  a  Taylor’s  series,  and  retaining  only  the  terms 
of  order  (AL)^  or  lower,  one  obtains  for  the  vectors  /f],  K2,  K^,  and  K^: 

K,  =  ALC(z)y(z)  (B-9) 

K2=ALC(z)y(z)  +  ^(mHC(z)fyiz)+j(ALfC'(z)y{z)  (B-10) 

=  (B-11) 

/fs  =  ALC(z)y(r)-H(AL)^(C(z))^y(z)  +  (AL)^  C'(z)y(z)  (B-12) 

Substituting  these  terms  into  Equation  B-8,  one  obtains: 

y(z -I- AL)  =  y(z)-i- ALC(z)y(z) 

-»-i(AL)^(C(z))^y(z)  (B-13) 

+ j(AL)^  C'(z)y(z) 
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The  next  step  is  to  substitute  the  explicit  expressions  for  y  and  C  (Equations  B-2  and 
B-3)  into  Equation  B-13.  This  yields; 


K„(z  +  AL)  =  V,(2)-;AZ,2,j3„/,(z)+AZ,J^r^V;(z) 


**1 


V.  i=i 


(B-14) 


The  term  Tm  (where  the  prime  denotes  differentiation  with  respect  to  z)  can  be 
evaluated  by  invoking  the  differentiation  theorem  or  surface  integ^  with  variable 
surface  (Ref.  20,  main  text): 


^  jJxdxdy-=  ^X(tan0)d/  (B-15) 

S(»)  5(1)  '  C(») 

where  ^  is  shown  in  Figure  1 1  and  C(z)  is  the  contour  defining  S(z).  The  relation  can 
be  used  to  show; 


=  ^^•e,(tan^)df  (B-16) 

S(«)  “  5(1)  C(i)“ 


Equation  B-14,  with  the  term  T'^  given  by  Equation  B-16,  gives  the  AL  and  (AL)^ 
terms  for  the  difference  equation  representation  of  the  Runge-Kutta  method.  The 
differences  between  this  method  and  the  boundary  discretization  technique  are 
embodied  in  terms  of  order  (AL)^  or  higher — the  AL  terms  are  exactly  the  same.  In  the 
limit  of  vanishing  AL,  the  differences  between  the  two  methods  vanish  and  both 
Equations  7  (from  the  main  text)  and  B-14  reduce  to  Equation  3  of  the  main  text. 
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APPENDIX  C. 

TRAVELING  WAVE  FORMULATION  OF  THE 
WAVEGUIDE  TRANSITION  PROBLEM 


The  forward-traveling  wave  coefficient  and  the  backward-traveling  wave 
coefficient  A^YzJ  are  related  to  the  equivalent  voltage  and  equivalent  current  of  the  pth 
mode  via  the  relations 


V^=^(A/  +  A-)  (C-1) 

and 

/p=Vi7(V-V)  (C-2) 

where  the  modal  wave  admittance  is  the  inverse  of  the  modal  wave  impedance  Z  . 
Using  the  boundary  discretization  approach,  the  equations  relating  the  wave  coeffi¬ 
cients  at  2  and  2  +  AL  (Fig.  1)  are; 


^  K  +  ( Z  +  AL)  = 


VW  (z)  +  ^Z^(z-hAL)  ^ 

r-i  (C-3) 


^^Yi{2  +  AL)  K^e--'''<‘"^>^'%-(z  + AL)  = 


t=i 


^^{2)e*'^'^'^A;(2)^'^^Yt{2+mV^ipe*>^’^'*^^'^A;{z  +  AL)  (C^) 


/=l 


where: 


=  JJef(z)  ep  dxdy 
s 


(C-5) 


In  implementing  Equations  C-3  and  C-4,  a  solution  is  obtained  and  scattering  matrices 
are  determined  according  to  Reference  7  (Chap.  5).  One  might  attempt  to  use  the 
equivalent  voltages  and  currents  in  conjunction  with  die  boundary  discretization 
approach  (Equation  7),  but  the  problem  with  this  scheme  is  that  the  complex  trigono¬ 
metric  functions  associated  with  the  evanescent  modes  will  quite  often  exceed  the 
numerical  range  of  the  computer.  In  this  study  this  problem  was  circumvented  by  using 
traveling  waves  according  to  the  procedure  described  by  Chu  and  Itoh.’ 
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The  differential  equations  for  the  wave  amplitudes  are  found  by  substituting 
Equations  C- 1  and  C-2  into  Equations  3  and  4.  These  equations  (from  Ref.  2  and  13) 
are; 


where: 


dA* 

p 

dz 


1  4(InZJ  _ 


dA 


p  _ 


dz 


^Pt 


'  2' 

dz 

1  d(lnZp)  ^ 

2 

dz  ' 

_  1 

r-1 


r-i 


(C-6) 


(C-7) 


(C-8) 


In  many  instances,  Equations  C-6  and  C-7  are  hard  to  implement  numerically 
because,  if  a  particular  mode  goes  through  cutoff  in  the  transition,  the  associated 
parameter  will  diverge  if  the  mode  is  TE,  and  the  associated  parameter  will 
diverge  if  the  mode  is  TM.  Therefore,  in  our  woik,  we  have  used  Equations  3  and  4 
instead  of  Equations  C-6  and  C-7. 
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APPENDIX  D. 

CALCULATION  OF  MEASUREMENT  ERRORS 


As  an  aid  in  establishing  rough  estimates  for  the  measurement  errors  in  our  tests, 
scattering  parameter  data  were  obtained  for  a  PRD  XIOIV  rotary  attenuator.  This 
attenuator  consists  of  a  rectangular  waveguide  with  a  movable  resistive  card  in  it  If 
this  card  is  oriented  perpendicular  to  the  modal  electric  field  lines,  little  or  no 
attenuation  should  occur.  Rotation  of  this  card  away  from  this  position  increases  the 
transmission  loss,  which  is  read  from  a  dial.  Figure  D- 1  shows  the  measured  S2 1  for  the 
settings  0.0  dB  and  0.2  dB.  The  difference  between  the  two  sets  of  measured  data  is 
0.2  dB,  ±  0.01  dB.  Therefore,  based  on  these  data,  we  make  the  claim  that  our  I52i  I  data 
are  accurate  to  within  dbO.Ol  dB. 

In  order  to  get  an  estimate  of  the  error  in  ISj  1 1,  we  performed  measurements  on  two 
rotary  attenuators  in  series.  First,  measurements  on  all  four  scattering  parameters  were 
made  for  both  phase  and  magnitude.  In  the  second  test,  a  short  was  placed  on  port  2, 
and  the  reflection  coefficient  was  measured.  We  denote  the  reflection  coefficient  in  the 
second  test5 1 1 .  This  parameter  iseasily  shown  to  be  related  to  the  scattering  parameters 
of  the  first  test  by: 


Su=Sn-^  (D-» 

I  +  O22 

The  relationship  between  the  actual  reflection  coefficients  and  the  measured  data  can 
be  written  as: 


5i'i 


(D-2) 


5jl  - 

(D-3) 

S22  =  6^22,^03 

(D-4) 

Since  all  three  measured  parameters  were  in  the  same  range  (30  to  35  dB),  it  seems 
reasonable  to  postulate  that  e  takes  on  the  same  value  in  all  three  relationships. 
Substituting  Equations  D-1,  D-2,  D-3,  and  D-4  into  Equation  D-1,  one  obtains  the 
following  equation  for  e: 

^22mtat(^U0ftar  ~  (^Umta/  ~  ^IZxuof  ^2liiuax  (D-5) 

Equation  D-S  has  two  roots,  one  of  which  oscillates  in  the  range  30  to  40  dB  between 
1 1  and  12  GHz.  Since  this  root  would  result  in  a  reflection  coefficients  near  0  dB  in 
Figures  S  and  7  (of  the  main  text),  it  can  be  rejected  as  unphysical.  The  other  root  is 
shown  as  a  function  of  frequency  in  Figure  D-2.  From  these  data,  the  calculated  error 
in  ISi  jl,  is  seen  to  be  within  ±10.0  dB.  Perhaps  not  coincidendy,  this  is  the  magnitude 
of  the  oscillations  in  the  data  of  Figures  S  and  7. 

Finally,  it  should  be  noted  that  these  error  bounds  are  somewhat  greater  than  those 
described  in  the  Wiltron  360  manual.  One  possible  reason  for  this  is  that  the  efficacy 
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of  the  calibration  procedure  can  be  degraded  due  to  faulty  cables  and  connectors  even 
if  the  network  analyzer  is  in  excellent  condition. 


Forward  Transmission 


Log.  Mag.  t>  Ref-o.OOOdB  0.100  dB/Div 


Figura  D-1 .  Measured  |Sj,|  for  a  precision  rotary  attenuator.  Top  line:  Attenuator  setting 
at  0.0  dB.  Bottom  line;  Attenuator  setting  at  0.2  dB. 


Figure  D-2.  Estimated  error  in  {SnI- 
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APPENDIX  E. 

NUMERICAL  ANALYSIS  OF  THE  STEEP  CIRCULAR 
WAVEGUIDE  TAPER 


Figure  E-1  shows  a  circular  waveguide  iris  with  steeply  tapered  sides.  Each  of 
these  two  sides  can  be  thought  of  as  a  linear  transition  of  length  L.  Experimental  results 
for  L  =  0  can  be  found  in  Ragan's  book  (see  the  Reference  list  following  the  main  text. 
Ref.  21,  pp.  213-214).  Figure  E-2  shows  these  experimental  results  along  with 
numerical  results  obtained  using  a  nine-mode,  mode-matching  routine.  The  two 
abscissas  in  Figure  E-2  represent  the  diameter  of  the  waveguide  aperture,  and  the  two 
ordinates  rq>resent  the  reflection  coefficient  5)t  magnitude  and  phase.  It  should  be 
mentioned  that  Carin  has  previously  investigated  this  device  and  has  published  results 
for  the  magnitude  of  5]  i  very  similar  to  our  own  (Ref.  8),  but  no  results  for  the  phase 
of  Si  1.  We  have  also  generated  results  for  L  =  OOlXo  and  L  -  llh^Xo  where  is  the 
free-space  wavelength  at  the  operating  frequency  9.375  GHz. 

These  data  were  generated  taking  into  account  nine  modes  and  using  two 
analytical  techniques.  First,  boundary  discretization  and  mode-matching  were  applied 
with  100  junctions  between  z  =  0  and  z-L.  Second,  the  Galeikin’s  method  version  of 
the  moment  technique^  was  used  with  six  cascaded  sections  and  five  triangle  weighting 
functions  per  section.  This  procedure  was  executed  at  six  iris  diameters  and,  at  each 
diameter,  the  results  for  5i  i  (mode-matching  vs.  method  of  moments,  L  -  0,00 1 X,,  vs. 
L  =  1O~^0)  were  indistinguishable;  these  data  are  shown  in  Figure  E-3  along  with  the 
L  =  0  mode-matching  results.  The  nonzero  L  results  converge  to  values  close  to  the  L 
=  0  results  in  magnitude  and  to  values  indistinguishable  from  theL = 0  results  in  phase. 


^  '  U- 


Figur*  E-1 .  Circular  waveguide  irisr  2d  •>  1 5/1 6  inch,  f  ■  1  /32  inch, 
frequency  «  9.375  GHz. 
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A. 


B. 


FIgur*  E-3.  Computsd  raftoction  coefficient  for  the  circular  waveguide  iris  of  Fig.  E-1 . 
The  solid  line  was  generated  for  the  L  -  0  case  using  mode-matching  with  nine  modes. 
In  addition,  data  were  generated  for  0.001  Xq  and  L-  10~*  Xq  i»ing  the  moment 
method  and  the  boundary  discretization  method.  The  same  results  (shown  using 
crosses)  were  obtained  for  bo<h  lengths  and  both  methods.  A.  Magnitude,  B.  Phase. 
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